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I ABSTRACT 

The  piston  modification  for  the  NAHBE  (Naval  Academy  Heat  Balance 

Engine)  consists  of  a cap  which  extends  above  the  standard  piston  used 

in  an  Otto  cycle. This  modification,  together  with  an  auxiliary  air 

inlet  between  the  intake  valve  and  fuel  supply,  has  been  shown  to  produce 

a significant  influence  on  combustion  within  the  cylinder.  One  theory 

to  explain  this  behavior  is  that  time  dependent  heat  addition  may  be 

affected  in  a field  of  pressure  waves  sustained  by  pressure  exchange 

(21 

between  the  two  chambers.'  ' 

The  objectives  of  this  study  were  twofold.  First,  the  transient 
convective  boundary  conditions  on  the  cap  were  studied  to  estimate  what 
probable  values  of  heat  transfer  coefficients  and  environmental  tempera- 
tures were  creating  the  temperature  behavior  within  the  cap  material. 
Second,  the  heat  transfer  to  and  from  the  cap  and  the  surrounding  gases 
within  the  cylinder  was  calculated. 

A mathematical  model  was  developed  which  predicts  the  transient, 
periodic  temperature  response  of  the  cap  as  a function  of  time  dependent 
convective  boundary  conditions.  Total  heat  transfer  to  and  from  the 
cap,  along  with  instantaneous  heat  flux  rates,  were  obtained  for  various 
operating  conditions.  Two  sets  of  convective  heat  transfer  coefficients 
were  used,  a set  of  low  values  and  a set  of  high  values.  The  actual 
convective  heat  transfer  process  is  not  fully  understood,  but  these 
initial  parametric  results  give  an  indication  of  the  regenerative  effect 
due  to  heat  transfer  from  the  cap  to  the  air-fuel  mixture  which  can  be 
expected  under  various  conditions. 
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II  FOREWORD 

The  aluminum  cap  on  the  piston  head  within  a NAHBE  engine  undergoes 
a complex  thermal  response  during  the  engine  operation.  The  complex 
surface-gas  interaction  during  combustion  is  not  fully  understood  and 
thus  the  resistance  to  heat  transfer  by  convection  and  radiation  between 
the  gas  and  solid  cap  can  only  be  estimated.  However,  parametric  studies 
using  ranges  of  value  which  include  the  suspected  actual  values  can  give 
a useful  insight  into  the  thermal  process. 

In  this  study  the  top,  bottom,  and  edge  surfaces  of  the  cap  are 
exposed  to  three  different  environmental  conditions  through  computer 
simulation.  The  three  environmental  temperatures  and  the  three  surface 
heat  transfer  coefficients  can  change  independently  during  each  stroke 
of  the  4-stroke  engine.  The  speed  of  the  engines  can  be  varied  and  re- 
sults are  given  in  Section  IV  for  speeds  between  500  and  5,000  rpm.  The 
simulation  model  gives  the  temperature  of  the  cap  as  a function  of  radius 
and  time,  the  instantaneous  surface  heat  transfer  rates,  and  the  total 
heat  transfer  to  and  from  the  cap  during  a specified  time  interval.  A 
listing  of  the  computer  programs  used  in  this  simulation  appears  in 
Appendix  A. 
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Ill  MATHEMATICAL  MODEL 


The  cap  is  modeled  as  a radial  fin  of  constant  thickness  (6).  The 
inside  base  radius  r^.  = 1.0  in.,  the  outside  radius  r^  = 1.75  in.,  and 
the  thickness  (initially)  6 = 0.01  in.  The  temperature  is  assumed  to 
be  uniform  across  the  fin  thickness  at  a given  radius  and  time.  The 
differential  equation  for  this  transient,  one-dimensional  heat  conduction 
problem  is  given  by: 


an 

3r^ 


ni?-  T«  - •'B  (T-T-.b)  - »£  = jIt 

Here,  hj,  hg,  and  hg  are  the  time  dependent  heat  transfer  coefficients 


on  the  top,  bottom,  and  edge  surfaces  respectively.  Likewise,  T 


,T’ 


Too,b,  ^00, E dependent  environmental  gas  temperatures 

for  the  three  respective  regions. 

The  boundary  conditions  are  as  follows: 


at  r = , T = T^  = 500  F 

at  r = r,.  - k(|I)^^  = h,  (T-T.^,) 

^0 


The  initial  condition  is: 

at  t = 0,  T = 500;  r^-  i r ^ r^  . 

An  explicit,  finite  difference  model  is  used  to  approximate  the 
above  partial  differential  equation.  A grid  consisting  of  thirteen 
nodes  (ar  = 1/16  in.)  is  chosen  to  approximate  the  fin.  Application 
of  the  first  law  of  thermodynamics  to  this  grid  gives  the  following 
set  of  explicit,  finite  difference  equations. 
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. 1 /aAtwArx  ,,  n ^ n, 
^ I ^aF-^^TT)  ^^-+1  - ^i-l’ 


* (^)  l''!  T.,t  " ''B  T.,bI 
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«xi  n hr-Ar 

T^-'’  = (1/G)  [T^._!f  + (-^) 

, ^ /hpAr/kv 

+ n - UE.^... ) 5 t; 


Here  represents  the  temperature  at  node  i (1  £ i _<  13)  for  a given 
time,  t = n.  T.  represents  the  temperature  at  node  i for  the  time 
t + At  = n + 1 . The  constant  G = (AK)^/aAt. 

Numerical  stability  requires  that  the  time  step  At  be  chosen 
such  that: 

<#>  ("b  * "t*  * 2(f^)  < 1.0 

Tor  an  aluminum  fin  with  h 'v  50  BTU/hr-ft^-  F,  this  '"eouires  At  < 0.01 

sec.  In  the  NAHBE  analysis  this  does  not  impose  a limitation  since 

much  smaller  time  increments  are  used  to  account  for  the  rapid  variations 

in  h and  during  the  4-stroke  cycle. 

The  assumed  time  variation  o^  T the  environmental  gas  temoera- 

«.T , 

ture  above  the  fin,  is  shown  in  Figure  1 for  five  complete  cycles  at  an 
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RPM  of  1500.  An  assumed  maximum  temperature  of  4500  F occurs  at  the 
end  of  the  expansion  stroke.  Intake  air  at  70  F and  exhaust  at  1218  F 
is  used.  The  gas  temperature  variation  beneath  the  cap  is  given  in  Figure 
2,  and  the  assumed  gas  temperature  variation  at  the  edge  is  shown  in 
Figure  3.  These  functions  can  be  easily  changed  as  a better  understanding 
of  the  actual  process  is  obtained. 

The  assumed  variation  of  the  convective  heat  transfer  coefficients 
with  time  are  estimated  based  upon  the  normal  range  of  expected  valves 
for  turbulent  flow  of  a gas  over  a surface.  In  order  to  establish  a 
benchmark  result  a low  value  of  h = 5 BTU/hr-ft^-F  was  used  for  hj,  hg, 
and  hg  during  intake.  A value  of  10  BTU/hr-ft^-F  was  used  during  the 
other  three  strokes  of  compression,  expansion,  and  exhaust.  Table  A 
gives  the  values  used. 


STROKE 

TOP 

EDGE 

BOTTOM 

EXHAUST 

10 

10 

10 

INTAKE 

5 

5 

5 

COMPRESSION 

10 

10 

10 

EXPANSION 

10 

10 

10 

Table  A:  Heat  Transfer  Coefficients 
for  Benchmark  Runs 


All  values  can  be  assigned  independently  as  more  knowledge  is  gained. 

The  benchmark  runs  all  used  the  above  values  of  the  convective  heat 
transfer  coefficients,  as  well  as  the  previously  described  environ* 
mental  temperature  variations. 

A second  set  of  computer  solutions  was  obtained  using  values  of 

convective  heat  transfer  coefficients  near  those  reported  in  Reference  3 

5 


OH  EDGE 


FIGURF.  3:  Variation  of  T , with  time 


based  upon  results  obtained  by  Eaton  Corporation  for  exhaust  valves. 
These  values  are  tabulated  in  Table  B. 


STROKE 

TOP 

EDGE 

BOTTOM 

EXHAUST 

60 

60 

40 

INTAKE 

80 

40 

30 

COMPRESSION 

60 

60 

40 

EXPANSION 

60 

60 

40 

Table  B: 

Heat  Transfer  Coefficients 

for 

Eaton  Runs 
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IV  DISCUSSION  OF  RESULTS 


A.  Benchmark  Series  with  Parametric  Variations 

Figure  4 shows  the  temperature  response  of  an  aluminum  cap  as 
a function  of  time  for  nodes  3,  5,  7,  9,  11,  and  13  using  benchmark 
conditions  for  h and  a speed  of  500  rpm.  The  initial  temperature  is 
500  F throughout  the  fin  and  a quasi-steady  state  condition  is  reached 
after  8 seconds.  As  can  be  seen,  node  3 reaches  a mean  temperature  of 
570  F with  oscillation  + 8°F.  Node  5 reaches  625  F,  node  7 reaches  665  F, 
and  node  9 reaches  685  F.  Nodes  11  and  13  are  affected  by  the  edge  con- 
vective boundary  condition  and  behave  in  approximately  the  same  manner. 

The  numerical  values  in  Table  1 show  that  the  edge  temperature  (node  13) 
oscillates  between  697  F and  705  F during  the  time  interval  8.7  l t £ 10. 
Table  1 also  shows  that  the  total  heat  transfer  from  the  cap  to  the  gas 
during  the  initial  10  second  interval  is  0.1093  BTU  and  the  total  heat 
transfer  from  the  gas  to  the  cap  in  1.73  BTU  during  the  same  interval. 

Figure  5 shows  the  transient  response  for  the  benchmark  run  at  1500  RPM, 
and  Table  2 gives  the  numerical  values  afte''  quasi-steady  conditions 
are  achieved.  The  mean  temperatures  obtained  at  1500  RPM  are  very  close 
to  those  at  500  RPM.  However,  the  temperature  oscillations  are  much 
smaller  at  1500  RPM  as  expected.  The  total  heat  transfer  from  the  cap  to 
the  gas  is  only  slightly  higher  at  1500  RPM.  e.g.  0.1118  vice  0.1093. 

The  5000  RPM  benchmark  results  are  given  in  Figure  6 and  Table  3. 

Note  the  scale  change  for  the  ordinate  in  figure  6.  Most  noticeable  are 
that  the  nodal  temperatures  have  increased  slightly  and  the  nature  of 

the  oscillations  have  radically  changed,  "he  maximum  temperature  reached 
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Table  3;  Numerical  Values  for  Transient  Response  at  5000  RPM 


''■4 

e''' 


at  the  edge  of  the  fin  is  25  F above  the  lower  RPM  valves.  The  heat 
transfer  to  the  gas  from  the  cap  is  now  0.1209  BTU  and  the  heat  transfer 
to  the  cap  is  1.822  BTU  during  the  first  ten  seconds  of  run  time.  When 
integration  is  carried  out  between  10  sec  ^ to  20  seconds,  with  the 
cap  temperatures  relatively  constant  over  the  entire  interval,  these 
heat  transfer  results  are  increased  by  approximately  20%  during  10  sec 
_<  t _<  2D  sec. 

To  show  the  parametric  effect  of  changing  the  convective  heat  trans- 
fer coefficient  during  intake,  a change  was  made  to  the  benchmark  values. 
Figure  7 shows  the  effect  of  increasing  the  convective  heat  transfer  co- 
efficients during  intake  on  the  top,  bottom,  and  edge  of  the  cap  by  a 
factor  of  5,  from  5 BTU/hr-ft^-F  to  25  BTU/hr-ft^-F.  All  other  variables 
remain  the  same  as  the  benchmark  runs.  At  500  RPM,  the  results  in 
Figure  7,  compared  with  Figure  4,  show  that  the  temperatures  along  the 
cap  are  reduced  by  approximately  40F.  Comparison  of  Table  4 with  Table  1 
shows  that  the  heat  transfer  from  the  cap  to  the  gas  is  now  0.4809  BTU 
rather  than  0.1093  BTU,  an  increase  by  a factor  of  4.4,  The  heat  transfer 
to  the  cap  is  also  slightly  higher,  1.775  BTU  vice  1.730  BTU.  It  mist 
be  remembered  that  these  results  all  depend  upon  the  somewhat  arbitrary 
assumption  of  the  gas  temperatures  shown  in  Figures  1-3. 

Figure  8 and  Table  5 show  the  same  effect  for  the  same  high  intake 
h values  at  1500  RPM.  Again,  a 40  F temperature  drop  is  observed  and 
the  same  increase  in  heat  transfer  to  the  gas  is  obtained. 

Figure  9 and  Table  6 are  included  to  show  that  approximately  the 
same  effect  can  be  obtained  by  increasing  only  the  convective  heat  transfer 
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FIGURl:  8:  Transient  Response  at  1500  RPM  for  lligli  h Intake  Values 
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TABLE  5:  Numerical  Values  for  Transient  Response  at  1500  RPM 
for  High  H INtake  Values 
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TABLE  6:  Numerical  Values  for  Transient  Response  at  1500  RPM  for 
High  h Intake  on  Top  Cap  Surface  Only 


coefficient  on  the  top  surface  of  the  cap  to  25  BTU/hr-ft^F,  and  leaving 
an  others  at  the  benchmark  values.  This  can  be  seen  by  comparing  Fig- 
ures 8-9,  and  Tables  5-6.  In  this  case  the  total  heat  transfer  from  the 
cap  to  the  gas  is  0.4363  BTU  during  the  first  ten  seconds. 

Instantaneous  heat  transfer  values  are  also  shown  in  the  last  column 
of  the  previous  six  tables.  When  the  base  cap  temperature  is  kept  at 
500  F,  as  assumed,  the  heat  transfer  is  predominantly  positive  (into 
the  cap  from  the  environment).  Figures  10-12  give  the  same  results  in 
graphical  form  for  three  runs  at  1500  RPM.  As  pointed  out  earlier,  the 
greatest  amount  of  negative  heat  transfer  (from  the  cap  to  the  gas) 
occurs  when  all  convective  heat  transfer  coefficient  values  during  intake 
are  increased  by  a factor  of  5.  This  is  shown  in  Figure  11. 

All  results  shown  in  Figures  1-12  are  for  an  aluminum  cap.  The 

properties  of  aluminum  used  for  the  simulation  were  3.3  ft^/hr  for  thermal 
diffusivity  and  155  BTU/hr-ft-F  for  thermal  conductivity.  For  comparison 
a benchmark  run  at  1500  RPM  was  made  for  a steel  cap.  Properties  of  the 
steel  were  0.48  ft^/hr  and  26  BTU/hr-ft-F  for  thermal  diffusivity  and 
thermal  conductivity  respectively. 

Figure  13  gives  the  transient  response  for  steel  at  1500  RPM  with 

benchmark  values.  Comparisons  with  an  aluminum  cap  at  1500  RPM  (Figure  5) 

lead  to  the  expected  conclusions.  The  thermal  oscillations  in  the  steel 

are  much  more  damped  than  in  aluminum,  and  the  time  constant  for  steel  is 

much  larger  than  for  aluminum.  In  fact,  after  13  seconds  the  edge  of  the 

cap  has  exceeded  1000  F and  is  still  increasing.  Thus,  the  use  of  a steel 

cap  increases  its  operating  temperature  by  more  than  300  F over  most  of 
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BENCHMARl  1500 


FIGURE  10:  Instantaneous  Meat  Transfer  at  1500  RPM 

with  Benchmark  values 


the  cap  surface,  as  compared  to  an  aluminum  cap.  The  heat  transfers 
to  the  gas  from  the  cap  during  the  first  ten  seconds  at  1500  RPM  is 
0.1712  BTU  for  a steel  cap  as  compared  to  0,118  BTU  for  an  aluminum  cap. 
This  increase  is  due  to  the  increased  operating  temperature  of  the  cap. 
Conversely,  the  heat  transfer  from  the  gas  to  the  cap  is  1.538  BTU  for 
the  steel  and  1.728  BTU  for  the  aluminum. 

B.  Eaton  Series  with  Parametric  Variations 

Results  for  computer  simulations  using  the  convective  heat 
transfer  values  given  in  TABLE  B begin  with  Figure  14.  Comparison  of 
Figures  14  and  4 show  that  the  fin  reaches  quasi-steady  temperatures 
faster,  and  the  predicted  temperatures  are  much  higher.  Aluminum  fin 
edge  temperatures  around  1,000  F are  predicted  by  the  model.  The  heat 
transfer  to  the  gas  from  the  fin  is  increased  from  0.109  BTU  to  2.63  BTU, 
a significant  change! 

At  1500  RPM  the  difference  between  the  benchmark  values  and  Easton 
values  are  shown  in  Figure  5 and  15.  For  this  case  the  fin-gas  heat 
transfer  increases  from  0.111  BTU  to  2.87  BTU.  Similar  trends  were  noted 
at  5000  RPM.  The  corresponding  instantaneous  figures  are  Figures  6 and 
16.  The  heat  flux  values  for  the  benchmark  runs  (Figures  10-12)  show 
much  less  negative  heat  transfer  (from  fin  to  gas)  than  in  the  Eaton  runs 
(Figures  17-19).  Of  course  the  positive  heat  transfer  is  also  greater. 

When  a steel  fin  is  used  with  the  Eaton  convective  heat  transfer 
coefficient  values,  temperatures  at  the  edge  of  the  fin  exceed  1200  F 
as  shown  in  Figure  20.  The  heat  transfer  from  the  steel  fin  to  the  gas 
during  the  first  ten  seconds  is  3.7  BTU. 
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TABLE  7:  Numerical  Values,  Eaton  Series,  RFM 
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Figure  18:  Instantaneous  Heat  Flux  at  1500  RPM 
with  Eaton  Valves 
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Instantaneous  Heat  Flux  at  5000  RPM 
with  Eaton  Valves 
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For  a parametric  change  to  the  Eaton  values  of  h,  all  values  of  h 
were  increased,  to  80  BTU/hr-ft^-F.  The  results  for  1500  RPM  are  shown 
in  Figure  21  for  an  aluminum  cap  and  should  be  compared  with  Figure  15. 

The  heat  transfer  to  the  gas  for  this  case  was  3.9  8TU  during  a ten 
second  interval.  For  a steel  cap  with  the  same  high  h values,  the  response 
is  shown  in  Figure  22  and  the  heat  transfer  obtained  was  5.4  BTU. 

To  show  the  effect  of  cap  thickness  the  thickness  was  doubled  while 
keeping  all  h values  at  80  BTU/hr-ft^-F.  The  effect  of  increasing  the 
cap  thickness  under  these  conditions  can  be  seen  by  comparing  Figures  23 
and  21.  The  maximum  temperature  in  the  cap  is  reduced  from  1096  F to 
958  F.  The  heat  transfer  to  the  gas  for  the  thicker  fin  is  3.0  BTU 
rather  than  3.8  BTU  during  the  first  ten  seconds. 
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Figure  22:  Transient  Response  of  a Steel  Cap  with  aVi  h valves 
set  at  80  BTU/hr-ft^-F 
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Figure  23:  Transient  Response  for  a Thick  Cap 
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V CONCLUSIONS 

The  mathematical  model  and  resulting  computer  simulation  appear  to 
give  a reasonable  prediction  of  how  the  NAHBE  engine  cap  behaves  from  a 

V heat  transfer  point  of  view.  Further  verification  await  measurements  on 

\ * 

the  actual  engine  during  operation.  Changes  to  either  the  assumed  values 
..  of  the  convective  heat  transfer  coefficients  or  the  time  dependent  func- 

j tions  which  describe  the  temperature  of  the  surrounding  gases  on  top, 

'\ 

bottom,  or  edge  of  the  cap  will  give  results  which  are  significantly 

V different  from  those  presented  in  this  report. 

II 

Observations  resulting  from  these  computer  simulations  are  as 
. follows: 

a.  The  benchmark  runs  give  temperature  distributions  within 
the  aluminum  cap  which  appear  reasonable,  although  no  experimental  data 
is  available  at  this  time.  However,  the  heat  transfer  from  the  cap  to 
the  gas  is  too  small  to  create  a significant  regenerative  effect. 

b.  The  Eaton  runs  give  heat  transfer  values  which  are  between 
10-20^  of  the  estimated  work  of  compression  ( m C^aT)  and  thus  can  produce 

a regenerative  effect.  However,  the  predicted  cap  temperatures  are  probably 
higher  than  what  actually  occurs. 

c.  Replacement  of  the  aluminum  cap  with  a steel  cap  increases 
the  cap  operating  temperature  by  approximately  300  F. 

d.  Increasing  the  thickness  of  an  aluminum  cap  by  a factor  of 
2 decreases  the  operating  temperature  by  100  - 150  F. 

' If  regenerative  heat  transfer  is  to  be  important  while  operating 

temperatures  are  reasonable,  the  model  suggests  that  the  heat  transfer 
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coefficient  (e,g,  non-dimensional  gas  temperature  gradient  at  the  cap 
surface)  must  be  low  when  the  environmental  gas  temperature  is  high  (to 
reduce  heat  transfer  to  the  cap),  and  high  when  the  environmental  gas 
temperature  is  lower  than  the  cap  temperature  (to  increase  heat  transfer 
from  the  cap).  One  hypothesis  is  that  during  compression  the  filling  of 
the  cavity  beneath  the  cap  prevents  the  build  up  of  a "thick"  boundary 
layer  on  the  top  cap  surface  and  this  produces  a much  higher  heat  transfer 
from  the  cap  to  the  gas  during  compression.  Also,  it  is  possible  that 
during  the  complicated  combustion  process,  gradients  may  be  set  up  in 
the  ionized  gases  to  effectively  reduce  the  heat  transfer  coefficient, 
or  the  effective  temperature  difference  between  surface  and  gas,  during 
the  combustion/expansion  stroke. 


43 


REFERENCES 


Blaser,  R. , Pouring,  A.,  Keating,  E.,  and  Rankin,  B.,  "The 
Naval  Academy  Heat  Balance  Engine,"  USNA  Report  EW  8-76, 

June,  1976. 

Pouring,  A,  A.,  et  al , "The  Influence  of  Combustion  with  Pressure 
Exchange  on  the  Performance  of  Heat  Balanced  Internal  Combustion 
Engines,"  paper  770120,  Society  of  Automotive  Engineers, 

February  28,  1977, 

Davis,  L.  J.  "Engine  Valve  Cooling  for  Emission  Control,"  Paper 
730055,  presented  at  The  Society  of  Automotive  Engineers, 

January  8-18,  1973. 


44 


I 


1 

I 


' I 

m 


VI 

APPENDIX  A 
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11:25:56 


PACE  1 


100  LIPWPyiNmAL*.“EllP’."l3UT‘ 

110  FILE  #1; ‘TRANS' 

120  FILE  #2; ‘HEAT* 

130  FILE  #3:  "INIT‘ 

140  SCRATCH  II 
150  SiCRATCH  #2 
160  DIM  0(3000) 

170  LET  01=1 
ISO  DIM  T(20).F(20) 

190  DIM  P(20) 

200  REM  INPUT  CAP  THICKNESS! IN) 

210  LET  D=0. 01 

220  REM  INPUT  CAP'  LEMJTHdN) 

230  LET  L=0.  75 

240  REM  INPUT  WD  OF  NODES 

25<)  LET  N=13 

260  PRINT' INPUT  DELTA  TIME"; 

270  INPUT  10 
230  LET  )(=L/(N-1) 

290  PRINT 

300  REM  INPUT  THERMAL  DIFFU3IVITY(FT^2/HR) 

310  LET  A=3. 3 

320  REM  irWT  THERMAL  C(»JI»JCTIVITY(BTU/}«-FT-F) 

330  LET  K=155 

340  LET  F=A*I0H44/((X^2)t36OO) 

350  PRINT 

360  PRINT'INPUT  PRINT  INTERVAL(SEC),  MAX  TIMEISEO'i 
370  INPUT  P,  E 
380  PRINT 

3?'0  REM  AS:;iCN  INITIAL  TEMP  DISTRIBUTION 
400  CALL*INITIAL":N,L,K.P().T(),i3 
410  LET  Pl=f 

420  REM  APLY  EXPLICIT  FINITE  DIFFERENCE  EONS 

430  CALL'EXP":  D,  N, L, A, K,  F, P,  E,  T( ),  0( ). 01, HI, H2,  H9, 10,  #2,  #3 

440  REM  PUT  TRANSIENT  RESPONSE  ON  FILE  #1 

450  LET  Z9=INT(N/2)+l 

460  LET  23=INT(E/P1) 

470  FOR  1=1  TO  19 

480  FOR  J=1  TO  Z8+1 

490  PRINT  #l:Ji','iO((J-l)*Z9+I> 

500  NEXT  J 

510  PRINT  #1;'1E37,1E37' 

520  NEXT  I 

530  PRINT  I2;“1E37.1E37' 

540  END 
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PAGE  1 


1 


r 

100 

SUB"  INIT1AL":N,LK.P(),T  0.13 

110 

«AT  T=ZER(N+1) 

120 

LET  X=0 

1 

130 

FOR  1=1  TO  13 

140 

’H'^UT  I3:P(I) 

XT  I 

. 

140 

KiR  1=1  TO  N 

170 

LET  T(I)=P(I) 

180 

LET  X=X+Ly(N-l) 

e 

190 

NEXT  I 

200 

SUBEND 

< 

r 


c 
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PAGE 
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9 


iw 

no 

i;o 

130 

140 

150 

160 

170 

ISO 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

■230 

340 

350 

360 

370 

380 

390 

40'3 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 


S'JD“EXP“:  D.  N,  L,  A.  K.  F,  P,  E,  T( ),  0( ).  01.  HI.  H2.  H9. 10.  #2.  #3 
REM  DEF  H ON  TOP 
['EF  FNX(ZO) 

IF  Z0<X1  THEN  180 
IF  Z0<X2  THEN  200 
IF  Z0<;X4  THEN  220 
LET  FNX=10 
GO  TO  230 
LET  FNX=25 
GO  TO  230 
LET  FNX=10 
GO  TO  230 
LET  FNX=10 
FNLND 

DEF  FNB(ZO) 

REM  DEF  T GAS  ON  BOTTOM 
IF  Z0<X1  THEN  350 
IF  20<X2  THEN  310 
IF  Z0<X4  THEN  330 

LET  FNEs( (W5-W4)/(X5-X4) )*(Z0-X4)+H4 
GO  TO  360 

LET  FNE=((H2-«1)/(X2-X1))»(Z0-X1)+W1 
GO  TO  360 

LET  FNB=((U4-W3)/(X4-X3))«(Z0-X3)+«3 
GO  TO  360 
LET  FNB=«0 
FNEND 

DEF  FNT(ZO) 

REM  DCF  T GAS  ON  TOP 
IF  ZOai  THEN  480 
IF  Z0<X2  THEN  440 
IF  Z0<X4  THEN  460 
LET  fNT=Y2 
GO  TO  490 

LET  FNT=( ( Y2-Y1 )/(X2-Xl 1 )*(Z0-X1 )+Yl 
K)  TO  490 

LET  FNT=((Y4-Y3)/(X4-X3))»(Z0-X3)+Y3 
GO  TO  490 
LET  FNT=Y0 
FfENO 

DEF  FNY(ZO) 

REM  DEF  H ON  BOHOH 
IF  Z0<X1  THEN  570 
IF  Z0<X2  THEN  590 
IF  Z0<X4  THEN  610 
LET  FNY=10 
GO  TO  620 
LET  rNY=5 
GO  TO  620 
LET  FW=IO 
GO  TO  620 
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EXP  0S/;:/77  11:27.17  PAGE  2 


« 610  LET  FNY=10 

620  FNEND 

»»■’  630  DEF  FHE(ZO) 

Vj  C 640  REM  DEF  T GAS  ON  EDGE 

65(»  ZOCXl  THEN  700 
7'  Z0<X2  THEN  720 

Z0<X4  THEN  740 
630  LET  FKE=1000 
690  GO  TO  750 

>,  C 700  LET  FNE=530 

> 710  GO  TO  750 

j.'  720  LET  FNE=300 

" C 730  GO  TO  750 

740  LET  FNE=1213 
f 750  FNEND 

‘ C , 760  DEF  FNZIZO) 

770  REM  DEF  H ON  EDGE 
'S')  IF  Z0<X1  THej  330 
4;-  7'?0  IF  Z0<X2  THEN  850 

* SijO  IF  Z0<X4  THEN  870 

810  LET  FNZ=10 
l,  C S-I'O  GO  TO  880 

' 830  LET  FNZ=5 

" TO  ^ 

' FNZ=10 
uU  TO  880 
370  LET  FNZ=10 
880  FNEND 
890  REM  INPUT  RPM 
900  LET  N5=500 
i 910  LET  X9=60/N5 

920  REM  DEF  X-Y  FOR  T TOP 
■930  LET  X0=0 
i 940  LET  Y0=70 

950  LET  X1=0.  5*X9 
960  LET  Yl=530 
C 970  LET  X2=X9 

980  LET  Y2=1218 
990  LET  X3=X2 
^ 1000  LET  Y3=4000 

1010  LCT  X4=3»X9/2 
1020  LET  Y4=4500 
C 1030  LET  X5=2#X9 

1040  LET  Y5=Y2 
1050  REM  DEF  Y FOR  T BOTTOM 
^ 1060  LET  H0=530 

1070  LET  Ul=530 
1080  LET  W2=1218 
^ l -'^O  LET  H3=1218 


no  LET  W4=1218 
1110  LET  U5^41 
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DIM  U(30).X(20),R(20) 

DIM  Q(20) 

REM  CALC  DELTA  R 

LET  X6=U(N-1)  ' : 

' REM  CALCllATE  INITIAL  VALUES  > ■ • 

LET  Z=0 
LET  Z0=0 
LET  H1=FNX(Z0) 

LET  H2=FNY(Z0) 

LET  H?=FNZ(ZO) 

LET  S9=(3. 1415?*(1.  75''2-l)/144)t(HU(FNT(IO)-T(l))+H2*<FNB(ZO>-T<l)J)+lH9*3, 14159*3. 25*0. 01/144)*(FNT(Z0)-Ta)) 

) CALL“0UT": N. L.  A. K,  Z.  T( ). 0( ). 01, S9, #2 
) REM  A°PLY  FINITE  DIFFERENCE  EONS 
) LET  C=P 

I LET  U5=A*I0*12/(K*D»3600) 

I LET  U(1)=T(1) 

I LET  B=H9*X5/(12*K) 

I LET  09=0 
I LET  P9=0 
I LET  N9=0 
) LET  H1=FNX(Z0) 

) LET  H2=FNYIZ0)  , - 

) LET  H9=FNZ<Z0) 

) LET  R<1)=1 
) FOR  1=2  TO  N-1 
) LET  R{I)=1+(I-1>*X6 

) LET  U(I)=T(U*(l-U5*(Hl+H2)-2*F)+F*(Tn+l)+T(I-l))+(F/2)*(X6/R(I»*(T(I*l)-T(I-l))+U5*(Hl«FNT(Z0)+H2«FNB(Z0)) 

> LET  Q(I)=(3. 14l59*(R(ir2-R(I-l)''2)/144)*(Hl*(FNT(Z0)-U(I))+H2*(FN6(Z0)-U(I))) 

) NEXT  I 
) LET  1=1+1 

) LET  U(I)=F*(T(I-l)+B*FNE<ZO))+(l-f*U+B))*T(I) 

) LET  Q(I)=(3. 14159*(1.  75''2-R(N-l)*2)/144)*(Hl*(FNT(Z0)-U(Il»+H2*(FNB(Z0)-U(IJ))+«9*(3. 14159*3.25*0  01/144)*(F>C1Z0)  U 
3 LET  S9=0 
3 FOR  1=2  TO  N 

3 LET  S9=S9+Q(I)  ’ 

3 NEXT  I 

3 IF  S9<0  THEN  1510 
3 LET  F9=P9+S9*I0/3600 
3 GO  TO  1520 
) LET  N9=N9+S9*I0/3600 
3 LET  09=Q9+S9*10/3600 
) LET  Z=Z+IO 
) LET  I5=INT(Z/(2*X9)) 

) LET  Z0=Z-I5*2*X9 

) IF  Z<C  THEN  1600  . . 

) CALL'OUT'.N, LA, K.Z, TO, 00,01,89, #2  J 

3 LET  C=C+P 
3 IF  Z>E  THEN  1640 
3 FOR  1=2  TO  N 
3 LET  T(I)=(J(I) 

3 lEXT  1 

50 
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^ 1630  GO  TO  1320 

1640  PRINT 

1650  PRINT'TOTAL  Q=“i 09; "AFTER";  Z; “SEC*; “ TINE  INC=“;  10 
T'  1660  PRINT 

1670  PRINT"NEG  Q=*;N9;*P0S  Q=*;P9 
1^^''  '^TCH  »3 
C ''-0  FOR  1=1  TO  13 

i.vO  PRINT  #3;U(I) 

1710  NEXT  I 
O 1720  SUBEND 


o 


OUT 
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O 100  SUB’0UT»:N.L.A,K.Z,T().O().Ol,S9.i2 

110  [‘IH  0(201. F(20hQ<20I.X(20) 
120'CALCULATE  AND  PRINT  TINE 
C 130  LET  T=Z 

140  LET  T=INT(T»l(r3+0.  5)/10^3 
150  PRINT  T, 

O 160'CALaiLATE  DIMENSIONAL  TEMPERATURE 

170  FOR  1=1  TO  N 
180  LET  D(I)=T(I) 

C I’O  LET  F(I)=INT(D(I)*10+0.5)/10 

200  NEXT  I 

210' PRINT  ODD  NODAL  TEMPERATURES 
C 220  FOR  1=1  TO  N STEP  2 

230  PRINT  F(I)i 
240  LET  0(01)=F(I) 

C 250  LET  01=01+1 

260  NEXT  I 
270  PRINT“HT=’iS9 
C 280  PRINT  #2:  Z.‘. ‘,89 

2?0'CALaiLATE  SURFACE  HEAT  FLUX 
300  LET  0(1)=K*(D(2)-D(1))/(L/IN-1)) 

G 310  LET  X(1)=0 

320'CALD.'LATE  INTERfW.  HEAT  aUXES 
330  FOR  1=2  TO  N-1 

-340  LET  0(I)=K*(D(I+l)-D(I-l))/(2*iy(ll-l)) 
350  LET  X(I)=X(I-l)+L/(N-l> 

360  NEXT  I 
370  LET  e(N)=0 
380  LET  X<N)=L 

390'PRINT  HEAT  aUX  VALl£S  IF  DESIRED 
O 400  REM  PRINT  LOOP 

410  SU6EN0 


